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ABSTRACT 

Context. Ice desorption affects the evolution of the gas-phase chemistry during the protostellar stage, and also determines 
the chemical composition of comets forming in circumstellar disks. From observations, most volatile species are found in 
\ H20-dominated ices. 

Aims. The aim of this study is first to experimentally determine how entrapment of volatiles in H2O ice depends on ice 
thickness, mixture ratio and heating rate, and second, to introduce an extended three-phase model (gas, ice surface and 
Oh ' ice mantle) to describe ice mixture desorption with a minimum number of free parameters. 

Methods. Thermal H20:C02 ice desorption is investigated in temperature programmed desorption experiments of thin 
(10 - 40 ML) ice mixtures under ultra-high vacuum conditions. Desorption is simultaneously monitored by mass spec- 
i-C ■ trometry and reflection-absorption infrared spectroscopy. The H20:C02 experiments are complemented with selected 
H20:CO, and H20:C02:CO experiments. The results are modeled with rate equations that connect the gas, ice surface 
JL I and ice mantle phases through surface desorption and mantle-surface diffusion. 

Results. The fraction of trapped CO2 increases with ice thickness (10-32 ML) and H20:C02 mixing ratio (5:1 - 10:1), but 
' ■ not with one order of magnitude different heating rates. The fraction of trapped CO2 is 44 - 84 % with respect to the initial 
CO2 content for the investigated experimental conditions. This is reproduced quantitatively by the extended three-phase 
model that is introduced here. The H20:CO and H20:C02:CO experiments are consistent with the H20:C02 desorp- 
tion trends, suggesting that the model can be used for other ice species found in the interstellar medium to significantly 
t-H ■ improve the parameterization of ice desorption. 
> ' 
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; 1. Introduction jPontoppidan et al.ll2008|) . Such an ice desorption scheme 
OO . provides most of the gas phase reactants for the chem- 
. In pre-stellar cores, cold outer protostellar envelopes and istry taking plac e at later stages in these warm regions 
t-h ; protoplanetary disk midplanes, most molecules, except for (jDoty et al.N 2004L It is therefore crucial to understand ice 
, H 2 , are frozen out on dust grains, forming ice mantles. m i xtu re desorption and to effectively implement it in as- 
■ The main ice component in most lines of sight is H 2 0, trochemical networks. The aim of this study is to provide 
£ ; followed by CO and C0 2 , with a typical abundance of a laboratory basis for this process and to demonstrate how 
>^ ■ (0.5 - 1.5) x 10~ 4 for H?Q ice with respect to H 2 around it can be modeled both in the laboratory and in space. 
; solar-type protostars {van Dishoec^|2006|). Infrared obser- Laboratory experiments have provided most of the 
vations of pre-stellar cores show that most CO? ice and current knowledge about ice thermal desorption, in- 
some of the CO ice is mixed with H 2 jKnez et al.ll2005|). eluding desorption energies fo r mos t pure simple ices 
The remaining CO and C0 2 are found in separate ice (Sandford & Alla mando"laTll988l 119901: jFraser et all [2 001: 
layers. Based on these observations, H 2 and C0 2 are i Collings et all 12004 jOberg et all l2005t iBrown & Bo lina 
thought to form simultaneously on the grain surface dur- 20071: iBurke & Brownir2010|) ~Desorption from ice mixtures 
ing the early stage of cloud formation. When the cloud differs from pure ice desorption because of different bind- 
becomes denser, gas phase CO freezes out on top of the i ng energies between the mixture components (e.g., the 
water-ric h ice, resulting in a bi-lav ered ice mantle, as de- CO binding energy increases from 830 K in pure ice 
scribed in|Pontoppidan et al.| (|2008|). to 11 80 K in H 2 0-dominated ice mixtures l|Collings et all 
Once the pre-stellar core starts collapsing into a pro- 2003)) and because of trappin g of volatile spec ies in the 
tostar, it heats its environment, including the icy grains. H 2 hydrogen-bonding ices (Collings et al. 2004). Volatile 
This results in the desorption of the CO-rich layer into the components therefore desorb from H 2 0-rich ice mixtures 
gas phase, in structural changes in the water-rich ice layer, at a minimum of two different temperatures, correspond- 
and eventually in the desorption of the water-rich layer ing to the desorption of the species from the surface of the 
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H2O ice and from molecules trapped inside the bulk of the 
H2O ice, which only start desorbing at the onset of H2O 
desorption. Additional desorption is sometimes observed 
at the temperature for pure volatile ice desorption and 
during ice re-structuring, e.g., at the H?Q ph ase change 
from amorphous to crystalline jViti et alj|2004)) . This H2O 
restructuring occurs at <~140 K in the laboratory (for as- 
trophysical timescales the re-structuring temperature and 
desorption temperatur e decrease), which is close to the on- 
set of H 2 desorption (jCollings et alJl2004|) . 

Of the different ice mixture desorption features, the en- 
trapment of volatile species in H2O ice is astrochemically 
the most important to quantify. The trapping of CO in a 
water ice results in a factor of five increase in the effec- 
tive desorption temperature. In a recent cloud core col- 
lapse model, this corresponds to trapped CO desorbing at 
30 AU from the protostar compared to pure CO ice des- 
orbing at 3000 AU. The case is less dramatic, but still sig- 
nificant, for CO2, which desorbs at ~300 AU when pure, 
and at 30 AU if trapped in H 2 ice jAikawa et alj|2008t 
IVisser et alj |2009). Efficient ice trapping may therefore al- 
low some volatiles to stay frozen on the dust grains during 
accretion o f envelope mater ial onto the forming protoplan- 
etary disk jVisser et alj|2009|) . 

There are only a few models th at have incorporated 
the effects of ice mixture desorption. ICollings et al. ( 2004) 
investigated the desorption of 16 astro physical l y rele vant 
species from H 2 Q :X 20:1 ice mixtures . I Viti et alJ d2004j) and 
IVisser et a l. (2009) used the results of Collin gs et al.l (t2004) 
to split up the abundance of volatiles in up to four differ- 
ent flavors, with different desorption temperatures. These 
correspond to the fraction of each ice desorbing at the pure 
ice desorption temperature, from a H2O surface, during 
H2O ice restructuring and with H2O, respectively. This 
approach has provided information on the potential im- 
portance of ice trapping for the chemical evolution dur- 
ing star formation. However, this model does not take 
into account specific ice characteristics such as ice thick- 
ness, volatile concentration and heating rate, on which the 
amou nt of trapped volatiles in t he water ice may also de- 
pend (Sandford & Allamandolal ll988|) . These characteris- 
tics need to be determined experimentally to correctly pa- 
rameterize step models, where such are sufficient to model 
ice desorption. Strong dependencies on e.g. ice thickness 
or concentration would however warrant the develop- 
ment of a more continuous parameterization of ice des- 
orption than the assignment of flavors. 

These dependencies are naturally included in a few 
ice mixture desorp t ion models of spec ific binary ices 
llCollings et al.l 12003b iBisschop et al.ll2006|) . The molecular 
specificity of these models, together with a large number 
of fitting parameters has, however, prevented their incor- 
poration into larger astrochemical models. Therefore, in 
most gas-grain networks, desorption is still treated as if 
ices were pure, disregarding y olatile entrapment in the 
water matrix (e.g. lAikawa et al]|2008() . 

Another problem with current gas-grain codes is that 
evaporation is often incorporated as a first-order pro- 
cess, while it is experimentally found to be a zeroth- 
order process with respect to the total ice abundance for 
ices thicker than one monolayer. Desorption models from 
the last decades ha ye shown the necessity of using a ze - 
roth order kinetics jFraser et al. 2001; Collings et al. 2003). 
Incorporating ice desorption as a first-order process with 



respect to the total ice abundance effectively means that 
molecules throughout the whole ice are allowed to desorb 
at the same time, wh ich is non-physical l lFraser et all 2001; 
Bisschop et al] |2006|) . This can be solved by treating the 
bulk and surface of the ice as separate phases as it has been 
done by IPontoppidan et all (|2003|) and [Pontoppida n et al.l 
( 2008) for CO ice desorption and by ICollings et alJ (|2005) 
for H2O desorption and crystallization. Its successful use 
in astrochemical models makes this approach an attractive 
option to parameterize laboratory ice desorption, since 
the results can then be easily transferred into an astro- 
physical context. In this family of models, molecules are 
only allowed to desorb from the surface, which is contin- 
uously replenished by molecules coming from the man- 
tle, and therefore the desorption kinetics are automati- 
cally treated correctly. This model results in a zeroth-order 
desorption behavior, in agreement with the experiments, 
since the number of molecules available for desorption re- 
mains constant in time. The model also results in trapping 
of volatiles in the bulk of the ice since the mantle molecules 
cannot desorb into the gas phase. Th e three-phase model 
we bu ild on was first introduced by Hasegawa & HerbstI 
(|1993 ), but despite its advantages in treating different ice 
processes, it has not been generally used for ice mixture 
desorption, nor has it been further developed, presumably 
because it did not correctly reproduce the experimentally 
observed amount of volatiles trapped in the water ice. 

The goals of the present study are first to experimen- 
tally characterize how the trapping efficiency of CO2 in 
H2O ice depends on different ice characteristics (with com- 
plementary experiments on CO and tertiary mixtures) and 
second to use these experiments as a guide to improve our 
understanding of the trapping process within the three- 
phase model framework. The description of the extended 
three-phase model is explained in Section [2] The exper- 
iments used to get information on the volatile entrap- 
ment and to calibrate the model are described in Section 
[3] Laboratory results on H20:C02 ices, complemented by 
H20:CO and H20:C02:CO ice desorption results, are pre- 
sented in Section [4] Section [5] presents the model fitting 
parameters and model results. Finally, the consequences 
of treating ice mixture desorption with the extended three- 
phase model under astrophysical conditions are discussed 
in Section[6] 

2. Desorption Model 

This study addresses the desorption of volatiles mixed 
with water and how to predict the fractions of volatiles 
in the ice and gas phase during a warm-up of the ice. 
The model is a sys t em o f rate equations based on the 
Hasegawa & Herbst (1993) model, but with the addition 
of diffusion. It aims at providing a solution for the amount 
of volatiles trapped in water with respect to the ice char- 
acteristics that can be directly inclu ded into astrochemica l 
models, as used by lViti et alJ | |2004|) and lVisser e~ al. (2009). 
The model applies to species in the water-rich ice layer; 
the interface with an upper CO-rich ice layer is not treated 
here. 

2.1. Basic three-phase model 

The model used here to predict the trapping of volatile 
species in a water dominated ice is based on the three- 
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Fig. 1. Cartoon defining the ice mantle (white), ice sur- 
face (gray), and g a s ph ase (black) according to the 
Hasegawa & Herbst ( 1993) three-phase model. Panels a) 
and b) show the different phases before and after a des- 
orption event. 



phase model by Hasegawa & HerbsJ l|1993|) . In this model, 
gas-grain interactions are addressed by considering three 
phases: the gas phase, the surface of the ice and the 
bulk/mantle of the ice (Fig. 1). The original model in- 
cludes reactions between species in both gas and solid 
phase, as well as accretion from the gas to the ice and 
thermal and non-thermal desorption. Thermal desorption 
alone is presented here. The model is based on the princi- 
ple that molecules can only desorb from the surface into 
the gas phase and that the mantle molecules can only mi- 
grate to the surface following the desorption of a surface 
molecule. The time-dependent gas abundance of species i 
is given by 



dnf 



R 



where, 



R 



evap 



ve 



evap 



-Ei/T 



(1) 



(2) 



with nf and n| the gas phase and surface abundance of 
species z, respectively, E, its binding energy in K, and v a 
pre-exponential factor taken equal to 1 12 s" 1 , which is a 
standard value for physisorbed species ( Biha m et al.l2 001). 
The surface abundance of species i can be written as: 



dnf 
~dt 
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evap 
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repl 



where, 



R 



repl 
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ve~ E i /T 



(3) 
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where nf is the mantle abundance of species i, Yl nf 1 the 
total number of molecules in the mantle and a is the ice 
coverage on the surface, which is set to 2 ML to account 
for surface roughness. The first term in Eq. [3] represents 
the loss of molecules i from the surface by thermal desorp- 
tion. The second term is related to the replenishment of the 
surface sites by mantle molecules: the empty sites created 
by the desorption of any type of species from the surface, 

Yj |ve" E ' / ' T J nj, are statistically filled by molecules com- 
ing from the mantle. The probability for these molecules 
to be species of type i is equal to its mantle fraction, ■ 



The mantle abundance, nf, of species i changes according 
to 



dn T : 



df 



R 



repl- 



(5) 



Because of the term j^nr, the replenishment of the 

surface phase by the mantle molecules during ice mix- 
ture desorption depends only on the mixing ratio of each 
species in this model, e.g., for a H20:C02 1:1 ice mixture, a 
molecule that desorbs into the gas phase has a 50% chance 
to be replaced by a water molecule and a 50% chance to 
be replaced by a CO2 molecule. This results in desorption 
of some volatile species around the pure ice desorption 
temperature and the rest remains trapped in the water ice 
since water molecules quickly saturate the surface phase. 

The ice abundances and nf are all in cm~ 3 , a unit di- 
rectly related to the gas phase abundance. The abundance 
of species i on the surface is defined via the relation 



Nfn c 



(6) 



where Nf is the average number of molecules i on the 

grain surface, and zz d is the dust abundance. The same re- 
lation applies for the mantle abundance. 

2.2. Extended three-phase model 

The original three-phase model does not account for the 
preferred replenishment of the surface phase by volatile 
mantle species or that volatile species may diffuse more 
easily in the ice compared to water. Obere ~et al.l (|2009b|) 
showed that this diffusion can result in segregation of 
the ice components, which is important for temperatures 
well below the desorption energy of most volatile species 
in an ice mixture. This demixing mechanism changes the 
surface replenishment pro babilities proposed in the origi- 
nal three-phase model by Hasegawa & Herbst (1993). Our 
proposed extension of the three-phase model accounts 
for this by introducing a mantle-surface diffusion term. 
Trapping of volatiles still occurs, but the surface-mantle 
diffusion of volatiles is enhanced compared to the original 
model, resulting in that more than 50 % of the empty sites 
are filled by volatiles species. Quantitatively, this changes 
the surface and mantle abundances n s - and nf 1 as follows: 



dn s t 
~dT 



-R, 



evap 



+ R. 
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for i 7^ H2O, and where e^O-/ * s me ener gy barrier for 
a volatile molecule i and a water molecule to swap ( i.e. 
change position) within the ice and fj a fraction between 
and 1 that is described below. The expression for the gas 
phase abundance remains unchanged (see Eq. [J}. The dif- 
fusion term Rf lSi is added to the surface abundance (sub- 
tracted from the mantle abundance) to enhance the man- 
tle to surface circulation for a volatile species i at the ex- 
pense of the water, thus, Rf lii is expressed differently for 
volatile and water molecules. The volatile ice diffusion 
rate depends on the balance of the probability of volatile 
molecules to move from the mantle to the surface at the 
expense of a water molecule and on the probability of 
the reverse process. This swapping process probability de- 
pends on the energy barrier q 5 ^ of the process and on 
the energy difference before and after the swap, equal to 
Eh 2 o — £/• The diffusion rate for the water molecules is 
the negative sum of the diffusion rates for the volatiles, 
^H*0 = ~~ ^» lff ' smce m e total abundance of molecules 
in the mantle and in the surface is not affected by the 
diffusion process. A similar formalism was used to de- 
scribe H20:CC>2 segregation in lOberg et alJ (|2009b|) : an ex- 
change of a surface H2O molecule and a mantle volatile 
is generally energetically favorable because H2O forms 
stronger bonds than volatile species and a mantle H2O 
molecule can form more bonds compared to a surface H2O 
molecule. 

From segregation studies of binary ices, it has become 
clear that only a limited fraction of the mantle participates 
in the mantle-surface circulation and th at this fraction de - 
pends on the initial ice mixture ratio llOberg et alJl20 09b). 
This is represented by the fraction /, 

/, = !_-! -^LlJ_ (11) 

where «™' ml is the number of mantle molecules i ini- 
tially in the ice, c, an empirical factor determined for each 
volatile i, and xf 11 the initial mixing ratio of volatiles i 
with respect to water. The expression c,(xj m )^ describes 
the number of mantle molecules available for segregation 
for a particular ice mixture before the onset of desorp- 
tion and follows the form found by Oberg et al. ( 2009Q) in 
ice segregation experiments when f> is set equal to 2. The 

term n™' im — c^x™ 1 )^ is the number of mantle molecules 
protected from segregation. When the later expression ex- 
ceeds the current number of volatile mantle molecules nf 1 , 
fi reaches zero and segregation stops, i.e., the diffusion of 
volatile mantle molecules to the surface stops. Thus this 
definition results in a gradual slowdown of the 'upward' 
mantle-surface diffusion of volatile species, regulating the 
trapping characteristics of H2O ice for different volatiles. 

We have tested the performance of this extended three- 
phase model on the desorption of mixed H20:C02 ices 
by comparing model and experimental TPD experiments, 
where the model TPDs are constructed using the rate 
equations [TJ and [8] In the model TPDs the initial ice tem- 
perature is raised in steps proportional to the heating rate 
and at each time step the rate equations from the three- 
phase model are applied to calculate the temperature de- 
pendent desorption and diffusion rates. The desorption 
rate of the volatile is what is plotted in the TPD curves. 



TPD experiments of pure ices are performed to deter- 
mine the binding energies E,-. The other free parameters 
that are used to optimize the model are the swapping en- 
ergies e^j between H2O molecules and volatiles i and 
the empirical factor c, used to parameterize the diffusion 
of volatiles i from the mantle to the surface. These two pa- 
rameters are determined by performing TPD experiments 
of binary ice mixtures of H20:CC>2 with different mixing 
ratios, thicknesses and heating rates and by comparing the 
output of the model with the experimental trends, i.e., the 
amount of volatile species that remains trapped in the wa- 
ter ice at temperatures higher than the desorption temper- 
ature of the volatile species. 

It is important to note that the model does not include 
the finite pumping speed during experiments. This will af- 
fect the derived desorption barriers and these are there- 
fore not meant to replace the ones derived from more de- 
tailed pure ice experiments in the literature. As long as the 
pumping rate is constant with temperature, excluding the 
pumping rate will not affect the determined ice fraction 
that desorbs at a low temperatures versus the fraction that 
desorbs with H2O. This is a reasonable assumption above 
the pure volatile ice desorption temperatures, where cry- 
opumping is no longer efficient. The derivation of c, and 
% o'-i should therefore not be affected by this simplifica- 
tion. 



3. Experiments 

3.1. Experimental parameters 

The experiments in this study are chosen to simultane- 
ously provide data directly relevant to ice desorption in 
different astrophysical environments (with different ices) 
and to construct a proof-of-concept model for ice mixture 
desorption. The focus is on CO2 desorption from H2O ice 
mixtures, one of the most important ice systems around 
protostars, with supporting experiments on CO desorp- 
tion. While interstellar ices are expected to be complex 
mixtures, it is still useful to investigate desorption from 
binary H20:volatile ice mixtures since the H20:volatile in- 
teractions are expected to dominate the desorption process 
in space, both because H2O is the major ice constituent 
and because H2O generally forms stronger bonds with it- 
self and with volatiles than volatiles do. This hypothesis 
has been further tested by performing TPD experiments 
of two tertiary H^OiCC^CO ice mixtures. 

The ice thickness and structure in the experiments are 
chosen to be as similar as possible to the existing obser- 
vational constraints on interstellar ices. Interstellar ices 
are estimated to be less than 100 monolayers (ML) thick 
from the maximum amount of oxygen available for ice 
formation. The experimentally grown ices are between 10 
and 40 ML, since it is only possible to quantify ice thick- 
nesses up to a certain limit (40 ML in our case) using 
reflec tion-absorption infrared spectroscopy l|Teolis et al.1 
2007). Information on ice structure in space is limited, but 
the lack of a water dangling vibration at 3700 cm -1 sug- 
gests a less porous ice than typically produced in the labo- 
ratory. We minimized the porosity of ice analogues by in- 
jecting gas pe rpendicularly to the cold surface when g row- 
ing the ices ((Stevenson et al.lll999tlKimmel et al.ll2001^ . 



Edith C. Fayolle et al.: laboratory H20:CC>2 ice desorption data 



5 



3.2. Experimental procedures 

All desorption experiments are performed with 
CRYOPAD. This set-up ha s been described in detail 
elsewhere dFuchs et al.l l2006|) . The set-up consists of an 
ultra high vacuum (UHV) chamber with a base pressure 
of ~ 10~ 10 mbar at room temperature. Ices are grown 
on a gold-coated substrate situated at the center of the 
chamber that can be cooled down to 16 K by a close cycle 
He cryostat. The relative sample temperature is controlled 
with a precision of 0.1 K using a resistive heating element 
and a temperature control unit. The absolute sample 
temperature is given with a 2K uncertainty. The system 
temperature is monitored with two thermocouples, one 
mounted on the substrate face, the other on the heater 
element. 

A fourier transform spectrometer is used for reflection- 
absorption infrared spectroscopy (FT-RAIRS) to record vi- 
brational absorption signatures of molecules condensed 
on the gold surface. The spectrometer covers 700 — 
4000 cm -1 with a typical resolution of 1 cm -1 and an av- 
eraged spectrum consists of a total of 256 scans. Ice evap- 
oration is induced by linear heating of the substrate (and 
ice) in TPD experiments. RAIR spectra are acquired simul- 
taneously to monitor the ice composition during the TPD. 
A quadrupole mass spectrometer (QMS) is positioned at 
4 cm, facing the ice sample to continuously analyze the 
gas-phase composition mass-selectively and to obtain des- 
orption curves of evaporating molecules during the TPD 
experiments. 

Mixtures and pure gas samples are prepared from 
13 C0 2 (Indugas, min 99% of 13 C), C0 2 (Praxair, 99% pu- 
rity), 13 CO (Cambridge Isotope Laboratories, 98% purity) 
and from gaseous water at the saturation pressure of a 
de-ionized liquid sample at room temperature. The de- 
ionized water is purified by three freeze-pump-thaw cy- 
cles. The samples are prepared separately, then injected in 
the chamber via an inlet pipe directed along the normal 
of the gold surface. In all gas samples, an isotopologue of 
CO was used to separate the QMS signal from background 
CO and N 2 . Similarly, an isotopologue of C0 2 was used 
to minimize the overlap in RAIR spectra between C0 2 ice 
and atmospheric C0 2 gas outside the UHV chamber. 

H 2 and C0 2 ice amounts are determined directly us- 
ing the RAIRS band strengths provided by Oberg lTt alJ 
(2009c, a) for CRYOPAD. From these measurements the ab- 
solute ice thicknesses are known within 50%. The relative 
ice abundance uncertainties are smaller, <~20%, and due to 
small band strength variations with ice composition and 
temperature. 

Table Q] lists the set of TPD experiments performed 
to calibrate and test the desorption model presented in 
Section [2] The TPD experiments begin with the deposi- 
tion of pure or mixed ice samples on the gold substrate 
cooled to 16 — 19 K, and continue with a slow heating of 
the ices at a constant specified rate until the desorption 
of the molecules from the surface is complete. The evapo- 
rated gas phase molecules are continuously monitored by 
the QMS. RAIR spectra of the ices are acquired before heat- 
ing to determine ice thicknesses and mixture ratios as de- 
scribed above. Spectra are also recorded during the warm- 
up as a second independent way to determine the ice com- 
position and to monitor eventual structure modifications. 



The infrared data are reduced by subtracting a local 
baseline around the molecular features. Mass spectromet- 
ric data are reduced by subtracting the ion current from 
species present in the background for each mass channel. 
Absolute yields cannot be directly obtained by the QMS 
since it is situated away from the ice sample (4 cm) and 
thus some of the desorbing molecules may get pumped 
away before detection. All QMS desorption rate curves 
are therefore normalized in such a way that the time- 
integrated desorption rate from the various species cor- 
responds to their infrared spectrally measured ice abun- 
dance at the beginning of each experiment. 



Table 1. Overview of the desorption experiments 



Exp. 


Sample 


Ratio 


Thick. 


Heat, rate 


Trapped C0 2 /CO ice % 








(ML) 


(K.min- 1 ) 


wrt. C0 2 /CO 


wrt. H 2 


1 


H 2 




24 


1 






2 


13 co 2 




6 


1 






3 


"CO 




6 


1 






4 


H 2 0:C0 2 


10:1 


12 


1 


62 


6.2 


5 


H 2 0:C0 2 


10:1 


19 


1 


75 


7.5 


6 


H 2 0:C0 2 


10:1 


32 


1 


84 


8.4 


7 


H 2 0: 13 C0 2 


5:1 


32 


1 


64 


12.8 


8 


H 2 0: 13 C0 2 


5:1 


18 


10 


62 


12.4 


9 


H 2 0: 13 C0 2 


5:1 


18 


1 


53 


10.6 


10 


H 2 0: 13 C0 2 


5:1 


10 


5 


44 


8.8 


11 


H 2 0: 13 C0 2 


5:1 


10 


1 


45 


9.0 


12 


H 2 0: 13 C0 2 


5:1 


10 


0.5 


44 


8.8 


13 


H 2 0: 13 CO 


10:1 


14 


1 


43 


4.3 


14 


H 2 0: 13 CO 


10:1 


25 


1 


47 


4.7 


15 


H 2 0: 13 CO 


5:1 


20 


1 


24 


4.8 


16 


H 2 0: 13 CO 


2:1 


13 


1 


9 


4.5 


17 
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4. Experimental analysis 

4.1. Complementarity of RAIRS and QMS 

Figure |2] illustrates the agreement between desorption 
curves derived from QMS and RAIRS data for C0 2 in 
a 5:1 water-dominated H 2 0:C0 2 ice, 18 ML thick and 
heated at 1 K.min -1 rate (Exp. 11). The upper left panel 
in Fig. |2] shows the C0 2 stretching band recorded at dif- 
ferent temperatures during warm-up: after ice deposition 
at 22 K, at 62 K where segregation is known to be efficient 
(Oberg et al. 2009b), during the first ice desorption peak 
around 79 K, in the temperature interval between pure 
C0 2 desorption and H 2 desorption, and during desorp- 
tion of the trapped C0 2 . The right panel shows the des- 
orption rate of C0 2 derived from the same experiment by 
mass spectrometry. The bottom panel presents the cumu- 
lative ice loss versus temperature for this experiment, ob- 
tained both by integrating the C0 2 mass signal with re- 
spect to the temperature, and by integrating the C0 2 in- 
frared signal recorded at specific temperatures. The er- 
ror bars on the infrared data are due to variable ice band 
strengths with temperature and composition. Within these 
uncertainties the fractional ice loss curves derived by in- 
frared and by mass spectrometry agree well; there seems 
to be only a small systematic offset for the 80 - 130 K range. 
This implies that the first RAIR spectrum of the ice after 
deposition can be used to derive quantitative results from 
the TPD experiments. 
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Fig. 2. The upper left panel presents the infrared CO2 
stretching features at specific temperatures during the 
warm-up of 18 ML of a H 2 0:C0 2 5:1 ice heated at 1 
K.min -1 . The right upper panel presents the desorption 
rate of CO2 for the same experiment obtained by mass 
spectrometry. The bottom panel shows the ice loss for this 
experiment obtained by infrared measurements (crosses) 
and mass spectrometry (solid line). 



Figure |2] also shows that there is evidence for some ice 
loss between the two main desorption peaks. The cumu- 
lative QMS and infrared spectroscopy signals match each 
other at these intermediate temperatures, which points to 
that the measurements trace actual ice desorption in be- 
tween the pure ice desorption event and the desorption 
of trapped volatiles. The implications of this ice desorp- 
tion process is discussed below, but it is important to note 
that this is not incorporated into the model framework 
and this may be a limitation to step-wise desorption mod- 
els, whether using our parameterization or any of the pre- 
viously published ones. Quantifying this process would 
require an additional set of experiments where the mass 
spectrometer is mounted closer to the substrate to allow 
for the detection of very low desorption rates. 

4.2. Desorption trends 

Figure [3] shows the desorption of CO2 from Fi20:CC>2 ice 
mixtures of different thicknesses (a), with different CO2 
concentrations (b), and heated at different rates (d). In ad- 
dition there are two CO TPD curves from Fi20:CO mix- 
tures with different CO concentrations (c). For reference, 
Fig. |3^) presents the TPD curves of pure CO, CO2, and 
H2O ice heated at 1 K min -1 . The fraction of trapped 
volatile is obtained by integrating the QMS signal for tem- 
peratures above 110 K and dividing it by the QMS signal 
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Fig. 3. Experimental CO and CO2 desorption curves a-d) 
during warm up of ice mixtures (offset for visibility) to- 
gether with pure CO, CO2 and H2O ice TPD curves e). The 
heating rate is 1 K.min -1 except for when specified other- 
wise in a), the total ice thickness and mixing ratio are listed 
in for each experiment. 



integrated over the entire 20-160 K range. The chosen tem- 
perature of 110 K is well below the onset of the second 
desorption peak and the volatiles that desorbed during 
the first CO2 or CO desorption peak are (almost) entirely 
pumped, though as discussed above there seems to be a 
low-level type of desorption occurring between the main 
desorption peaks. Whether due to finite pumping or actual 
desorption this results in a 10-20% uncertainty in the de- 
termination of the trapped fraction, i.e., the choice of tem- 
perature integration limits affects the estimated amount of 
trapped ice by < 20%. The trapped percentage of volatiles 
in each experiment, defined with respect to the initial 
volatile ice content, is reported in the second last column 
of Table [T]The last column of Table [1] presents the trapped 
abundance of volatiles species with respect to the initial 
H2O abundance. This value is less variable compared to 
the trapped amount of CO/CO2 with respect to the initial 
CO/CO2 abundance presented in the preceding column. 
Both Tableland Fig.|3]show that for CO2 and CO the per- 
centage of trapped volatile species in the H2O ice is highly 
dependent on the experimental conditions; the CO2 trap- 
ping fraction varies between 44 and 84% with respect to 
the initial volatile content. In the following subsections, we 
report and discuss these dependencies. 
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4.2.1. Thickness dependency 

Figure |2K presents the desorption of a H20:C02 10:1 ice 
mixture for different initial ice thicknesses and shows that 
the amount of trapped CO2 (desorption around 140 K) in- 
creases with ice thickness. In contrast the amount of CO2 
desorbing around 70 K is independent of ice thickness in 
the experimentally investigated range. This implies that 
only CO2 molecules from the top part of the ice are avail- 
able for desorption at the CO2 desorption temperature. 
This can be explained by either a highly porous ice that 
allows CO2 to "freely" desorb from the top layers or by 
diffusion from the top layers of the mantle phase to the 
surface. In both cases the surface is eventually totally satu- 
rated by water molecules, trapping the rest of the volatiles 
in the ice mantle. 



4.2.2. Mixing ratio dependency 

The dependence of the volatile trapping with the mixture 
ratio is presented in Fig. [3j)) for H20:C02 and in Fig. |3j:) 
for H20:CO ice mixtures. In both cases the trapped frac- 
tion decreases as the volatile to H2O ratio increases. In 
other words, the amount of pores exposed to the surface 
or the diffusion length scale of volatiles in the ice must 
increase with increasing volatile concentration. A similar 
dependency was noted in lOberg et alJ (|2009b|) when mea- 
suring segregation in ices. Increased diffusion may either 
be due to a gradually looser binding environment in the 
volatile-rich ices or a break-down of H2O ice structure 
in the presence of higher concentrations of volatiles. The 
thick H20:CO ice experiments (Exps. 14, 15 and 17) with 
initial mixture ratios of 10:1, 5:1 and 1:1 show a continu- 
ous decline of the trapping fraction, which suggests that 
either the ice becomes continuously more porous or that 
the diffusion path length increases gradually with volatile 
concentration. 



4.2.3. Molecular dependency 



Similarly to ISandford & A llamandoTal l|1990|) and 
Collings et al.| (|2004j) we find that the trapping effi- 
ciencies of the investigated CO2 and CO in H2O ice are 
radically different. When comparing Fig. [3J)) and c), it 
appears that CO is much more mobile than CO2 in the 
H2O ice as demonstrated by the higher trapping fraction 
of CO? compared to CO in s imilar H2O ice mixtures. 
Sandford & Allamandola (1990) explained this difference 
from a combination of different binding energies of CO 
and CO2 in H2O ice due to molecular size, shape and 
electronic differences. These binding energies may equally 
affect the probability of escaping through an ice pore or 
diffusing through the bulk of the ice. 

4.2.4. Heating rate dependency 

In Fig. |3ji), the heating rate of the ice is var- 
ied for a H 2 0:C0 2 5:1 ice of 10 - 12 ML between 
0.5, 1 and5K.min _1 . This does not appreciably affect the 
trapping efficiency of CO2 in the H2O ice and implies that 
the process responsible for exchanges between the sur- 
face and the mantle is fast compared to the experimental 
warm-up time. If this was not the case, a lower heating rate 
would have resulted in a smaller amount of trapped CO2, 
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Fig. 4. a: Desorption rate of CO2 (solid line) and CO 
(dashed line) from two tertiary water dominated ice mix- 
tures (Exps. 18 and 19) with a 1 K.min -1 heating rate. - 
b: Implemented three-phase model desorption rate for the 
same eperiments. 



since a slower heating means more time for the migration 
of mantle molecules to the surface. 

The lack of a heating rate dependency on the trapped 
amount of volatiles (same trapped percentage in Exp. 10, 
11 and 12) also implies that there is a rather sharp bound- 
ary between the molecules in the upper layers that can dif- 
fuse to the surface (whether through pores or bulk diffu- 
sion) and molecules deeper in the ice that cannot. Even 
if the volatiles deep in the ice can diffuse within the ice 
mantle, diffusion 'upwards' must quickly become slow as 
the surface layers saturate with H2O molecules or alter- 
natively all accessible pores have been emptied. This ex- 
plains that the amount of desorbing CO2 molecules at low 
temperatures is thickness independent (the H2O 'ice cap' 
will become impenetrable after a certain amount of CO2 
molecules have desorbed) and that entrapment efficiencies 
are unaffected by lower heating rates. 

4.3. Tertiary mixtures 

The desorption rates for the tertiary mixtures (Exps. 18 and 
19) are presented in Fig.|4^). Similarly to the binary experi- 
ments, the trapping of volatiles is more efficient for a lower 
volatile to H2O ratio. The desorption curves for CO2 are 
not affected by the presence of CO. Thus the CO-CO2 in- 
teraction does not have a significant impact on the amount 
of CO2 trapped within the water matrix. Overall the TPD 
curves resemble the addition of desorption curves from 
two separate binary mixtures, except for a small fraction 
of CO that desorbs with C0 2 in the H 2 0:C0 2 :CO 11:4:1 
mixture (lower panel in Fig. |4). It is unclear whether this 
desorption is due to a co-desorption of CO with CO2 or to 
a release of CO that has been trapped under a barrier of 
CO2 surface molecules. The observed similarity supports 
the use of binary ice mixtures as templates to study dif- 
fusion and desorption even though they are not directly 
representative of interstellar ice mixtures. 

When comparing the desorption of CO from a tertiary 
and a binary ice mixture with the same H20:CO ratio and 



8 



Edith C. Fayolle et al.: laboratory H20:C02 ice desorption data 



ice thickness, it appears that less CO is trapped in the ter- 
tiary mixture. This may be due to ice structure changes, 
as discussed above, or to shielding of CO from the sticky 
water molecules by CO2 molecules (the CO-CO2 bond 
is weaker than the H2O-CO one), lowering the CO dif- 
fusion barrier. The CO2 and CO desorption curves from 
the dilute tertiary mixture both contain a small additional 
peak around 150 K, only seen elsewhere in the thickest 
H20:C02 10:1 ice experiment. A similar do uble peak was 
noted in the 20:1 desorption experiments of Collines et alJ 
(2004). 



4.4. Ice diffusion mechanisms: pore versus bulk diffusion 

The main mechanism behind diffusion in the ice mantle 
is not known and may differ between different ices. Most 
previous studies have focused on diffusion in cracks and 
pores and pore collapse has been introduced to explain ice 
trapping. The observation that both CO and CO2 become 
trapped even though they partly desorb at their, very dif- 
ferent, pure ice desorption temperatures is difficult to rec- 
oncile with pore collapse as the main trapping mechanism, 
however. That is, it would imply efficient H2O pore col- 
lapse both at ~30 K and ~70 K. 

Even if pore collapse does not provide a complete ex- 
planation of why ices become trapped, some kind of in- 
ternal surface hopping may explain why molecules can 
diffuse out of the ice. In this scenario, the CO2 desorp- 
tion ice thickness dependency is due to that the pores and 
cracks that are open to the surface only go down to a cer- 
tain depth, in this case a few ML for Fi20:C02 5:1. The 
different CO and CO2 trapping efficiencies may then be ei- 
ther due to different ice structures or to CO desorbing eas- 
ier through pores compared to CO2. Pore diffusion may 
thus be consistent with these particular experimental re- 
sults, but ice diffusion is present also in othe r ices that are 
know n to be quite compact, e.g. CO ice ( Bisschop et al. 
2006). While it is possible that diffusion occurs through 
completely different mechanisms in different ices, the con- 
cept of ice bulk diffusion has the prospect of approximat- 
ing mixing and de-mixing processes in all kinds of ices, 
regardless of structure. 

In the bulk diffusion scenario, an amorphous ice 
is viewed like a very viscous liquid, whose viscosity 
decreases with the volatility of the ice molecules. In 
lOberg et alJ (200%) this was modeled as molecules swap- 
ping places with a barrier significantly higher than sur- 
face hopping. Volatile molecules will tend to swap their 
way towards the surface because it is energetically favor- 
able to have the molecules that form weaker bonds in the 
surface layer (where fewer bonds can be made due to the 
ice-vacuum interface). Trapping is explained by that when 
volatile molecules diffuse from the top mantle layers to the 
surface and desorb, the top mantle layers become satu- 
rated with H2O and therefore viscous enough to be im- 
penetrable. The low desoption rate of volatiles between 
the volatile and H2O ice desorption peaks would however 
suggest that under some experimental conditions, small 
amounts of volatiles can escape through this H2O barrier. 
More experiments are required to test under which condi- 
tions this is a reasonable approximation. The model pre- 
sented below is an attempt to include the most important 
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Fig. 5. x 2 contour plot for fitting the model parameters 
% O-CO (parameterizing the ability of CO2 to switch po- 
sition with a H2O molecule) and Cqq 2 (parameterizing the 
CO2 ice thickness where swapping is fast compared to the 
investigated heating rates) using the experimentally deter- 
mined amount of CO2 trapped in a binary H20:C02 ice 
(Exps. 4 to 12). These two parameters regulate the distribu- 
tion of the volatile molecules in the gas, surface and mantle 
according to Eq. [Til and [TOl of Section [2] 

features of this concept while still keeping the number of 
parameters low. 

5. Model parametrization and performance 

5.1. Parametrization 

The pure ice desorption energies are derived from fitting 
the three-phase model to the experimental pure ice des- 
orption curves with the results: Eh 2 o = 4400 K, Eqq 2 = 
2440 K and E m = 1010 K . The H ?Q value is lower than the 
one found in Fraser et al. (2001). This discrepancy is prob- 
ably due to a combination of that we use a single experi- 
ment, do not include the pumping speed and fix the pre- 
exponential f actor to v = lO 12 ^ " 1 for every species here. 
The value in iFraser et alJ (2001) should thus still be used 
when modeling the absolute desorption temperature. For 
the purpose of parameterizing the desorption fractions we 
prefer our value for the sake of consistency. The remain- 
ing model parameters £^ 2 ^_j and c, are obtained sepa- 
rately for CO2 and CO. This is done through a x 2 anal- 
ysis, where trapping fractions from the model are com- 
pared to those from the binary experiments for a grid of 
£ B^-i an< ^ c i values. The minimum x 2 value for CO2 is ob- 
tained for e^ a ^ H O =2500 K, c CO2 =20.5 ML, but the swap- 
ping energy for H20:C02, which is linked to the ability of 
CO2 to swap with H2O within the solid phase, is not well- 
constrained between 1600-3200 K (Fig. [5)- In contrast, the 
parameter cqo 2 related to available amount of CO2 from 
the ice mantle that can migrate to the surface, is well con- 
strained which suggests that even at laboratory time scales 
mixture desorption is mainly governed by how large a 
part of the mantle is eligible for swapping with the sur- 
face, rather than the swapping barriers. 
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Fig. 6. Simulated CO2 and CO TPD curves from H20:CC>2 
and H20:CO mixtures for different thicknesses, ratios, and 
heating rates a-d). Panel e) presents the simulated desorp- 
tion of pure H2O, CO2 and CO ice. This figure connects the 
model outputs to the experiments shown in Fig. [3] 

The CO experiments can be fitted with £qq-y[ 2 q = 960 
K and cqq= 80 ML, but these are based on only a few ex- 
periments and the inequalities ^ho^co < e H20-C0 2 anc ^ 
c CO > C C0 2 are a l° ne well constrained. 

5.2. Model performance 

5.2.1. Desorption trends modeling 

Figure [6] shows the simulations of the binary mixture des- 
orption using the optimized model parameters from the 
previous section. Generally the qualitative agreement is 
good and the model captures the trends that were ob- 
served experimentally. The exact shapes of the modeled 
and experimental desorption curves differ for several rea- 
sons. First, the model does not take into account the range 
of environments from which the molecules desorb, e.g., 
the H2O/ volatile fraction changes during the desorption 
process and even a pure ice has a range of different bind- 
ing sites. This affects both the position and the width of 
the peak. Second, the model does not consider the dif- 
ferent water ice structures present at different tempera- 
tures; water ice crystallizes around 140 K (see the des- 
orption peaks for water in Fig. |3^), which may affect the 
shape and position of the second desorption peak of the 
volatiles. Finally, the model does not include finite pump- 



ing speeds, which results in abundance tails in the exper- 
imental curves. While these effects may all be important 
under special circumstances the aim of the extended three- 
phase model is not to reproduce the experimental results 
perfectly. Rather, the goal is to capture the main character- 
istics of ice mixture desorption. 

The increase of the trapped amount of volatiles with ice 
thickness and the experimental observation that the same 
amount of molecules desorbs around the volatile pure des- 
orption temperature, regardless of the ice thickness, are 
reproduced in the model because diffusion between the 
mantle and surface is only allowed from a fraction of the 
mantle, fj (Eq.[TT]in Section|2]l which depends on the kind 
of volatile and mixing ratio with water. This fraction is in- 
dependent of the ice thickness. Thus the same amount of 
volatile molecules migrates to the surface regardless of ice 
thickness, followed by saturation of the surface phase by 
water molecules. The rest of the volatiles is trapped in the 
mantle until H2O desorption, thus the trapped fraction de- 
pends on the ice thickness. 

The observed concentration effect on the trapping ef- 
ficiency is reproduced by the model because the fraction 
of volatiles migrating to the surface depends on the mix- 
ing ratio of the volatile with respect to water, xf 11 . The 
lower the concentration of volatiles in the ice is, the smaller 
the fraction of volatile molecules make it to the surface 
and the more become trapped. The higher mobility of CO 
compared to CO2 is also reproduced by the model as the 
molecular paramater Cqq 2 is lower than Cqq. Thus more 
volatiles are able to diffuse to the surface in the case of 
H20:CO mixtures. 

Experimentally, a low desorption rate is sometimes ob- 
served between the volatile and H2O desorption temper- 
ature. In the model, the diffusion barrier energy is low 
enough that the diffusion process is complete before des- 
orption of the volatile takes place and therefore there is no 
desorption between the pure and H2O desorption peaks. 
Such a low diffusion energy barrier is needed to repro- 
duce that the trapping efficiency is insensitive to the heat- 
ing rate (within the explored heating rate range). There is 
probably a second diffusion process at play at these inter- 
mediate temperatures, which cannot be reproduced by the 
current, simple parameterization. 



5.2.2. Quantitative agreement 

Figure[7]compares the volatile trapping fractions obtained 
by the optimized H20:C02 model to those found experi- 
mentally. The error bars include the uncertainties due to 
the choice of the temperature from which we integrate the 
second QMS peak and are between 10-20% for the differ- 
ent experimental data points. The error bars on the model 
results originate from the uncertainties in the input ice 
thicknesses, mixing ratios and binding energies and these 
were obtained by varying the model input values within 
the experimental uncertainty ranges and then comparing 
the model results. In general, the uncertainty in the mix- 
ing ratio has the largest effect, resulting in model predic- 
tion uncertainties of ~15%. We conclude that CO2 desorp- 
tion from binary mixtures is quantitatively described by 
the model. 
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Fig. 7. Percentage of CO2 experimentally trapped vs the 
simulated one. The triangles are the trapped percentages 
obtained from the binary mixture experiments. The dia- 
monds are the model predictions for the two tertiary mix- 
tures. The one-to-one ideal relation is plotted as a dashed 
line. 
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Fig. 8. The amount of CO2 ice during ice warm-up at 1 K 
per 100 years according to the three-phase model, assum- 
ing two different initial H20:C02 5:1 ice mixture thick- 
nesses (left panel) and two different 20 ML ice mixing ra- 
tios (right pannel). 



5.2.3. Predictive power 

Figure 0J) shows the output of the model for the 
H20:C02:CO tertiary ice mixture experiments (Exp. 18 
and 19). These experiments were not used to constrain the 
model and are as a test of its predictive power. The con- 
centration dependency in these experiments is reproduced 
by the model; an increase in the concentration of volatiles 
leads to a decrease in the trapping fraction. 

For the higher concentration mixture, H20:C02:CO = 
11:4:1 (Fig. |4Jd, bottom panel), the model gives a CO des- 
orption peak around 70 K (corresponding to the pure 
CO2 desorption temperature), peak that was also experi- 
mentally observed (Fig. [4^, bottom panel). In the model, 
this peak results from the formation of free desorption 
sites on the surface due to desorption of surface CO2. CO 
molecules that are mixed with water migrate to and des- 
orb from the surface easily since the swapping and bind- 
ing energies are very low compared to the CO2 values. 

In addition to reproducing these qualitative trends for 
tertiary mixtures, the three-phase model also treats cor- 
rectly the desorption order of both CO2 desorption peaks. 
Finally Fig. [7| shows that the model provides a reasonable 
quantitative agreement between the predicted and exper- 
imentally determined amounts of CO2 trapped in the ter- 
tiary ice mixtures (black diamonds). This is very promis- 
ing for extending this proof-of-concept model to more 
species and more complex mixtures. 

6. Astrophysical implications 

Trapping of volatiles in H2O ice is a crucial parameter 
when predicting t he chemical evo lution during star and 
planet formation (jViti et al.l [2004). The modified three- 
phase desorption model is used here to test the effect of 
different initial ice compositions and ice thicknesses on 
ice mixture desorption. Ultimately the three-phase model 
should, however, be integrated in a protostellar collapse 
model to simulate the ice desorption accurately during 
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Fig. 9. Amount of CO2 ice during ice warm-up for a 
H 2 0:C0 2 5:1 ice 20 ML thick simulated by different 
models : the implemented three-phase model described 
here (black solid line), the o riginal three-phase model 
byjHasegawa & Herbstl (|1993) (green dashed line), the 
IViti et al.l (|2004|) model (blue dash-dotted line). A heating 
rate of 1 K per century is used in the two first models and 
desorpti on around a 5 solar masses protostar is presented 
from the lViti et aH (|2004)) model case. 



star formation. The prime advantag e of the three-phase 
mode l, as initially introduced by Haseeawa & Fterbst 
(|1993|) and extended here, is that it can treat surface and 
ice chemistry correctly, since it differentiates between sur- 
face molecules that can react with gas phase molecules, 
and mantle molecules that are protected from further pro- 
cessing. 

Figure [8] shows the amount of CO2 ice with respect to 
the original H2O ice abundance as a function of tempera- 
ture for different ice thicknesses and mixing ratios heated 
at 1 K per 100 years , typical for infall of material during 
protostar formation (Torgensen et al. 2005|). 
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The percentage of CO2 entrapment in a diluted water 
ice is significantly affected by the initial ice thickness and 
mixing ratio; 50% of the initial CO2 abundance is trapped 
in a 10 ML ice and 95% in a 100 ML ice. A similarly dra- 
matic difference is seen when assuming different initial ice 
mixtures: 64% of the CO2 stays trapped in the 5:1 ice and 
the fraction increases to 84% for the 20:1 ice. 

The treatment of these trapping dependencies is one 
of the key strengths of the extended three-phase des- 
orption model presented here. Figure [9] compares CO2 
ice desorption from a H20:C02 ice using the extended 
three-phase model, the o riginal thr ee-phase mode l by 
lHasegawa & Herbstl (|1993|) , and the IViti et alj (|2004|) as- 
trochemical network. Assuming a 20 ML thick H20:CC>2 
5:1 ice heated at 1 K per century, our model predicts that 
64 % of the initial CO ? will be trapped by the w ater ice, 
while the model by Hasegawa & Herbst ( 1993) predicts 
an 80 % trapping amount. This differen ce originates from 
the la ck of mantle-surface diffusion in Haseg awa & Herbstl 
( 1993). Its implementation is clearly important to correctly 
treat trapping of volatiles and to account for segrega- 
tion observed around protostars (Ehrenfreund et al. 1998; 
iPontoppidan et al.ll2008t ). 

All the CO2 molecules are predicted to be trapped 
by the water ice when simulati ng H20:CC>2 ice des- 
orptio n with the IViti et al.l j2004|) model. The IViti et al.l 
(2004) model assumes a different heating rate compared 
to the one used for the two three-phase models, but 
this only affects the desorption temperatures and does 
not affect the volatiles trapping fractions. Instead, the 
high trapping fraction is due to the fact that the model 
was parametrized based on a desorption experiment per- 
formed for a H20:CC>2 ice with a ratio of 20:1, which 
differs from the 5:1 - 4:1 ratio found in dense molec- 
ular clouds and protost ellar envelopes jKnez et alj 12005; 
IPontoppidan et al.ll2008|) . In the case of a H 2 0:C0 2 20:1 
ices, our model outp uts agree w ell with 100 % trapping 
fraction used by lVitiet al. (2004), since we find that more 
than 95 % of the CO2 is trapped by the water ice for 10 - 
100 ML thick ices. 

These different model predictions demonstrate the 
need for systematic laboratory studies when modelling ice 
desorption, since ice properties, such as ice thickness and 
mixing ratio, affect the desorption process. Even when us- 
ing desorption step functions, the size of the step cannot be 
accurately decided from a single experiment. Rather the 
investment of multiple experiments are needed, together 
with their efficient parameterization, to obtain versatile 
models of ice desorption for arbitrary initial conditions. 
Already for binary ice mixtures, this results in large exper- 
imental data sets. It is therefore reassuring that using bi- 
nary mixtures as templates for more complex ice mixtures 
results in approximately the correct trapping predictions. 



7. Conclusions 

Desorption from H20-rich ice mixtures is complex in that 
the amount of trapped ice depends not only on the species 
involved, but also on the mixture ratio and the ice thick- 
ness ; there is no constant fraction of volatile species 
trapped in a H2O ice. This complex behavior can be re- 
pr oduced by extending, the three-phase model introduced 
by lHasegawa & Herbst) <|1993|) . 



Using the H20:C02 ice system as a case study, we 
showed that a three-phase model that includes mantle- 
surface diffusion can reproduce the amount of trapped ice 
quantitatively in a range of binary ice mixtures. The ap- 
propriate input parameters for the H20:C02 system are a 
swapping energy ej^^'_ C Q 2 =2250 K and a molecular pa- 
rameter cco 2 = 20-5 ML, which describes from which ice 
depth diffusion to the surface can occur. 

In the model, the different CO2 and CO behavior can 
only be reproduced if ej^o-CO < £ H 2 0-C0 2 and c co > 
C C0 2 - This suggests that diffusion/ molecule swapping in 
H20-rich ices depends equally on the breaking of the H2O- 
volatile bond (parametrized here by the swapping barrier 
height) and on the mass /volume of the diffusing volatile 
(parameterized here by the ice thickness that the molecule 
can diffuse through). The experimental trends found for 
H20:CO and H20:C02:CO ice mixture desorption are con- 
sistent with the H20:C02 trends, which suggests that the 
three-phase model is generally appropriate to model ice 
mixture desorption. 

However the ice desorption process is implemented in 
astrochemical models, this study demonstrates that it is 
vital to understand how ice mixture desorption depends 
on the ice characteristics. The extended three-phase model 
naturally treats ice desorption with the right kinetic order 
and reproduces volatile entrapment. Its use in astrochemi- 
cal networks for grain-gas interactions should improve the 
predictions of gas-phase and grain-surface species abun- 
dances in astrophysical environments. 
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